Alternative strategies for deciphering the genetic architecture of childhood Pre-B acute lymphoblastic leukemia

La leucémie lymphoblastique aigüe (LLA) est une maladie génétique complexe. Malgré que cette maladie hématologique soit le cancer pédiatrique le plus fréquent, ses causes demeurent inconnues. Des études antérieures ont démontrées que le risque à la LLA chez l’enfant pourrait être influencé par des gènes agissant dans le métabolisme des xénobiotiques, dans le maintient de l’intégrité génomique et dans la réponse au stress oxydatif, ainsi que par des facteurs environnementaux. Au cours de mes études doctorales, j’ai tenté de disséquer davantage les bases génétiques de la LLA de l’enfant en postulant que la susceptibilité à cette maladie serait modulée, au moins en partie, par des variants génétiques agissant dans deux voies biologiques fondamentales : le point de contrôle G1/S du cycle cellulaire et la réparation des cassures double-brin de l’ADN. En utilisant une approche unique reposant sur l’analyse d’une cohorte cas-contrôles jumelée à une cohorte de trios enfants-parents, j’ai effectué une étude d’association de type gènes/voies biologiques candidats. Ainsi, j’ai évaluer le rôle de variants provenant de la séquence promotrice de 12 gènes du cycle cellulaire et de 7 gènes de la voie de réparation de l’ADN, dans la susceptibilité à la LLA. De tels polymorphismes dans la région promotrice (pSNPs) pourraient perturber la liaison de facteurs de transcription et mener à des différences dans les niveaux d’expression des gènes pouvant influencer le risque à la maladie. En combinant différentes méthodes analytiques, j’ai évalué le rôle de différents mécanismes génétiques dans le développement de la LLA chez l’enfant. J’ai tout d’abord étudié les associations avec gènes/variants indépendants, et des essaies fonctionnels ont été effectués afin d’évaluer l’impact des pSNPs sur la liaison de facteurs de transcription et l’activité promotrice allèle-spécifique. Ces analyses ont mené à quatre publications. Il est peu probable que ces gènes de susceptibilité agissent seuls; j’ai donc utilisé une approche intégrative afin d’explorer la possibilité que plusieurs variants d’une même voie biologique ou de voies connexes puissent moduler le risque de la maladie; ces travaux ont été soumis pour publication. En outre, le développement précoce de la LLA, voir même in utero, suggère que les parents, et plus particulièrement la mère, pourraient jouer un rôle important dans le développement de cette maladie chez l’enfant. Dans une étude par simulations, j’ai évalué la performance des méthodes d’analyse existantes de détecter des effets fœto-maternels sous un design hybride trios/cas-contrôles. J’ai également investigué l’impact des effets génétiques agissant via la mère sur la susceptibilité à la LLA. Cette étude, récemment publiée, fût la première à démontrer que le risque de la leucémie chez l’enfant peut être modulé par le génotype de sa mère. En conclusions, mes études doctorales ont permis d’identifier des nouveaux gènes de susceptibilité pour la LLA pédiatrique et de mettre en évidence le rôle du cycle cellulaire et de la voie de la réparation de l’ADN dans la leucémogenèse. À terme, ces travaux permettront de mieux comprendre les bases génétiques de la LLA, et conduiront au développement d’outils cliniques qui amélioreront la détection, le diagnostique et le traitement de la leucémie chez l’enfant.

Caractérisation et étude d'un élément régulateur du gène codant pour le récepteur à la vasopressine de type 2

Thèse réalisée en cotutelle avec l'Université Pierre et Marie Curie, Paris VI, France

Étude du rôle régulateur de la lamine B1 dans l’activation plaquettaire : base moléculaire de la thromboprotection chez les patients porteurs d'anticoagulant lupique et d'anti-lamine B1

Les anticorps anti-phospholipides (aPL), tels que les anticoagulants lupiques (LAC), sont associés au développement récurrent de thromboses chez les patients atteints du lupus érythémateux disséminé (LED). Il a été observé que des titres élevés d’auto-anticorps antilamine B1 (anti-LB1), chez des patients porteurs de LAC, diminuent le risque de ces manifestations thrombotiques. Toutefois, la relation existant entre la lamine B1 (LB1), les anti-LB1 et la thromboprotection n’est toujours pas expliquée. Dans cette étude, nous avons donc cherché à comprendre comment la LB1 et les anti-LB1 induisent cette thromboprotection. Nous avons testé les effets d'anti-LB1 purifiés et de LB1 recombinante sur l'activation des cellules endothéliales et des plaquettes. Nous avons été en mesure de déterminer que la LB1, contrairement aux anti-LB1, possède une activité anti-plaquettaire. En effet, la LB1 réduit l’activation et l’agrégation plaquettaires in vitro et in vivo. Cette activité est due à une liaison directe de la LB1 aux plaquettes, suivie par une internalisation rapide dans des vésicules de clathrine. Par co-immunoprécipitation, nous avons découvert que la LB1 interagit avec le récepteur de l’insuline situé sur la membrane plaquettaire. La liaison de la LB1 à ce récepteur entraîne vraisemblablement son internalisation et l'inhibition d'une des cascades de signalisation normalement induite par le récepteur de l’insuline, menant éventuellement à l’inhibition des fonctions plaquettaires. L’ajout d’anti-LB1 purifiés dans nos expériences a permis d'augmenter de façon significative la persistance de la LB1 dans les plaquettes, une observation confirmée par la détection de LB1 uniquement dans les lysats de plaquettes prélevées chez des patients anti-LB1 positifs. iv Nos résultats suggèrent que la LB1 prend part aux mécanismes régulateurs des processus d’hémostase chez des sujets sains et que la présence d’anti-LB1, chez les patients lupiques, prolonge la persistance de cet auto-antigène dans les plaquettes, les empêchant ainsi de s’activer. Ce mécanisme expliquerait la diminution du risque de thrombose chez les patients LAC positifs porteurs d’anti-LB1 circulants.

Étude du polymorphisme intra- et inter-spécifique du gène β-tubuline chez des espèces de champignons mycorhiziens à arbuscules en vue de développer des marqueurs moléculaires

Les champignons mycorhiziens à arbuscules (CMA), classés dans le phylum Glomeromycota, ne peuvent pas être facilement identifiés par la morphologie de leurs spores et leurs mycélia à l'intérieur ou à l'extérieur des racines de leurs hôtes. Ce problème fondamental d'identification rend l'étude de leur diversité, en particulier dans leur habitat naturel (sol et racine) extrêmement difficile. Les gènes ribosomaux ont été largement utilisés pour développer des amorces spécifiques et en inférer des arbres phylogénétiques. Cependant, ces gènes sont très polymorphes et existent en plusieurs copies dans le génome des CMA, ce qui complique l’interprétation des résultats. Dans notre étude, nous avons étudié le polymorphisme intra- et inter-spécifique du gène β-tubuline, présent en faible nombre de copies dans le génome des CMA, afin d’obtenir de nouvelles séquences nucléotidiques pour développer des marqueurs moléculaires. Les gènes β-tubuline amplifiés à partir de l'ADN génomique de cinq espèces du genre Glomus ont été clonés et séquencés. L’analyse des séquences indique un polymorphisme intraspécifique chez trois espèces de CMA. Deux séquences paralogues très variables ont été nouvellement identifiées chez les G. aggregatum, G. fasciculatum et G. cerebriforme. Aucun polymorphisme n’a été détecté chez les G. clarum et G. etunicatum. Toutes les séquences montrent la présence de deux introns hautement variables. La majorité des substitutions ont été localisées dans les exons et sont synonymes à 90%. La conservation des acides aminés suggère un niveau élevé de sélection négative sur le gène β-tubuline et nous permet de confirmer que les CMA représentent un ancien groupe fongique (400 million d’années). L’analyse phylogénétique, réalisée avec vingt et une séquences nucléotidiques du gène β-tubuline, a révélé que les séquences des Glomaceae forment un groupe monophylétique bien supporté, avec les Acaulosporaceae et Gigasporaceae comme groupe frère. Les séquences paralogues nouvellement identifiées chez les G. aggregatum et G. fasciculatum n'ont pas été monophylétiques au sein de chaque espèce. Les oligonucléotides ont été choisis sur la base des régions variables et conservées du gène β-tubuline. Le test PCR des amorces β-Tub.cerb.F/ β-Tub.cerb.R a révélé des bandes spécifiques de 401 pb pour les séquences paralogues du G. cerebriforme. Deux paires d’amorces ont été développées afin d’identifier les séquences du groupe nommé Tub.1. Les tests PCR nous ont permis d’identifier certaines séquences du groupe Tub.1. Une paire d’amorce β-Tub.2.F/ β-Tub.2.R nous a permis d’identifier certaines séquences paralogues du groupe nommé Tub.2. L’analyse d’autres gènes combinée à celle du gène β-tubuline permettra le développement de marqueurs moléculaires plus spécifiques pour l’identification de CMA.

L’impact du polymorphisme du gène E2 sur la quantification de la charge virale du VPH-16 dans les maladies précancéreuses du col utérin

Le VPH-16 de même que certains VPH, dont le VPH-18, causent le cancer du col utérin. Son intégration dans le génome humain pourrait être un marqueur de progression de l’infection. Les charges virales totale et intégrée sont présentement mesurées en quantifiant par PCR en temps réel les gènes E6 (RT-E6) et E2 (RT-E2-1) du VPH-16. Nous avons évalué l’impact du polymorphisme du gène E2 sur la quantification de l’ADN du VPH-16 dans des spécimens cliniques. Dans un premier temps, le gène E2 de 135 isolats de VPH-16 (123 appartenaient au clade Européen et 12 à des clades non- Européens) fut séquencé. Ensuite, un test de PCR en temps réel ciblant les séquences conservées dans E2 (RT-E2-2) fut développé et optimisé. Cent trente-neuf spécimens (lavages cervicaux et vaginaux) provenant de 74 participantes (58 séropositives pour le VIH, 16 séronégatives pour le VIH) ont été étudiés avec les trois tests E2 (RT-E2-2), E6 (RT-E6) et E2 (RT-E2-1). Les ratios de la quantité d’ADN de VPH-16 mesuré avec RT-E2-2 et RT-E2-1 dans les isolats Européens (médiane, 1.02; intervalle, 0.64-1.80) et Africains 1 (médiane, 0.80; intervalle, 0.53-1.09) sont similaires (P=0.08). Par contre, les ratios mesurés avec les isolats Africains 2 (médiane, 3.23; intervalle, 1.92-3.49) ou Asiatique- Américains (médiane, 3.78; intervalle, 1.47-37) sont nettement supérieurs à ceux obtenus avec les isolats Européens (P<0.02 pour chaque comparaison). Les distributions des quantités de E2 contenues dans les 139 échantillons mesurées avec RT-E2-2 (médiane, 6150) et RT-E2-1 (médiane, 8960) étaient statistiquement différentes (P<0.0001). Nous avons observé que les charges virales totale (odds ratio (OR) OR, 2.16 95% intervalle de confiance (IC) 1.11-4.19), et épisomale du VPH-16 (OR, 2.14 95% IC 1.09-4.19), mais pas la présence de formes intégrées (OR, 3.72 95% IC 1.03-13.4), sont associées aux néoplasies intraepitheliales cervicales de haut grade (CIN-2,3), et ce, en contrôlant pour des facteurs confondants tels que l’âge, le taux de CD4 sanguin, l’infection au VIH, et le polymorphisme de VPH-16. La proportion des échantillons ayant un ratio E6/E2 > 2 pour les femmes sans lésion intraépithéliale (7 de 35) est similaire à celle des femmes avec CIN-2,3 (5 de 11, p=0.24) ou avec CIN- 1 (4 de14, P=0.65). Le polymorphisme du gène E2 est un facteur qui influence la quantification des charges intégrées de VPH-16.

Étude du polymorphisme du gène majeur d’histocompatibilité de classe IIb (MHIIb) chez l’omble de fontaine (Salvelinus fontinalis)

Les molécules classiques du complexe majeur d’histocompatibilité de classe II (CMHII) sont des glycoprotéines de surface spécialisées dans la présentation de peptides, principalement dérivés de pathogènes extracellulaires, aux récepteurs des lymphocytes T CD4+ afin d’initier la réponse immunitaire adaptative. Elles sont encodées, avec celles du CMH de classe I, par les gènes les plus polymorphiques identifiés jusqu’à maintenant, avec plusieurs loci et une grande diversité allélique à chacun d’eux. De plus, le polymorphisme des gènes du CMHII n’est pas limité qu’aux séquences codantes. Il est également observé dans les promoteurs où on a démontré ses effets sur le niveau d’expression des gènes. La variation de la régulation d’un gène est considérée comme un facteur important et pour laquelle des modifications morphologiques, physiologiques et comportementales sont observées chez tous les organismes. Des séquences d’ADN répétées impliquées dans cette régulation ont été identifiées dans les régions non-codantes des génomes. D’un autre côté, la sélection par les pathogènes permettrait l’évolution et le maintien du polymorphisme des gènes du CMH chez les vertébrés. À ce sujet, plusieurs études ont montré l’implication de différents allèles du CMH dans la résistance ou la susceptibilité aux maladies. Cette étude avait pour objectifs de caractériser le polymorphisme du gène MHIIb chez l’omble de fontaine (Salvelinus fontinalis) et de documenter ses effets au niveau de la survie conférée par des allèles et/ou génotypes particuliers lors d’une infection, ainsi que sur la variation du niveau d’expression du gène dans différentes conditions. Dans une première partie, nous avons identifié un total de 6 allèles du gène MHIIb, désignés Safo-DAB*0101 à Safo-DAB*0601, qui montrent une grande similarité avec les séquences codantes provenant de poissons téléostéens et de l’humain. L’analyse des séquences du domaine b1 a permis de détecter l’effet d’une pression sélective positive pour maintenir le polymorphisme dans cette région de la molécule. Quatre de ces allèles ont été testés lors d’une expérience d’infection avec le pathogène Aeromonas salmonicida afin d’évaluer l’effet qu’ils pouvaient avoir sur la survie des poissons. Nous avons trouvé que l’allèle DAB*0101 était significativement associé à la résistance à la furonculose. En plus d’avoir été identifié chez les individus homozygotes pour cet allèle, l’effet a également été remarqué au niveau de la survie les poissons de génotype DAB*0101/*0201. À l’opposé, les facteurs de risque élevé obtenus pour les génotypes DAB*0201/*0301 et DAB*0301/*0401 suggèrent plutôt une association à la susceptibilité. Étant donné la faible fréquence à laquelle l’allèle DAB*0101 a été retrouvé dans la population, le modèle de la sélection dépendante de la fréquence pourrait expliquer l’avantage conféré par ce dernier et souligne l’importance de ce mécanisme pour le maintien du polymorphisme du gène MHIIb chez l’omble de fontaine. Dans une seconde partie, nous avons rapporté la présence d’un minisatellite polymorphique formé d’un motif de 32 nucléotides dans le second intron du gène MHIIb, et pour lequel un nombre exclusif de répétitions du motif a été associé à chaque allèle (69, 27, 20, 40, 19 et 25 répétitions pour les allèles DAB*0101 à DAB*0601 respectivement). L’expression relative de quatre allèles a été évaluée dans des poissons hétérozygotes aux températures de 6 ºC et 18 ºC. Les résultats indiquent que les allèles possédant un long minisatellite montrent une réduction de l’expression du gène d’un facteur 1,67 à 2,56 par rapport aux allèles qui en contiennent un court. De même, des allèles qui incluent des minisatellites de tailles similaires n’affichent pas de différence significative au niveau de l’abondance du transcrit aux deux températures. De plus, l’effet répressif associé aux longs minisatellites est amplifié à la température de 18 ºC dans des poissons de trois génotypes différents. Nous avons finalement observé une augmentation significative par un facteur 2,08 de l’expression totale du gène MHIIb à la température de 6 ºC. Ces résultats appuient l’implication des séquences d’ADN répétées dans la régulation de l’activité transcriptionnelle d’un gène et suggèrent qu’un minisatellite sensible aux différences de températures pourrait être soumis aux forces sélectives et jouer un rôle important dans l’expression de gènes et l’évolution des organismes poïkilothermes.

Caractérisation de Cks1, régulateur du cycle cellulaire, dans le cancer épithélial de l'ovaire

Le cancer épithélial de l’ovaire est le cancer gynécologique le plus létal. La survie à 5 ans est de 30-40% chez les patientes atteintes d’une tumeur invasive(TOV), comparativement à 95% chez les patientes diagnostiquées pour une tumeur à faible potentiel de malignité (LMP). Au laboratoire, l’analyse de l’expression des gènes de la micropuce à ADN HuFL d’Affymetrix a révélé que Cks1 est un gène dont l’expression varie entre les tumeurs LMP et TOV. En effet, ce régulateur du cycle cellulaire est surexprimé dans les tumeurs TOV par rapport aux tumeurs LMP. Nous avons donc déplété Cks1 dans des lignées cellulaires tumorales invasives du cancer de l’ovaire dérivées au laboratoire, soit la TOV112D et la TOV1946, en utilisant des shRNAs sous le contrôle d’un répresseur inductible à la tétracycline. Puis, nous avons dérivé des clones stables inductibles à la tétracycline. Les résultats obtenus nous indiquent que la déplétion de Cks1 n’a pas d’effet sur la prolifération et la migration cellulaires, ni sur la formation de structures tridimensionnelles in vitro. Ainsi, nous pouvons conclure que Cks1 ne joue pas un rôle clé dans la progression tumorale par rapport aux paramètres testés. Or, des études supplémentaires seraient nécessaires pour expliquer les différences biologiques observées entre les deux types de tumeurs étudiées, et justifier cette variation observée de l’expression de Cks1.

Polymorphisme de la capside des papillomavirus appartenant à l’espèce 9 des alphapapillomavirus

Le gène L1 encode pour la protéine majeure de la capside des papillomavirus humains (VPH). L’information relative au polymorphisme de L1 pour les types autres que VPH- 16 est jusqu’ici limitée. Cet ouvrage explore le polymorphisme de L1 en comparant les séquences des types phylogénétiquement apparentés VPH-31, -33, -35 à VPH-16. Des spécimens génitaux recueillis de 732 femmes VIH-séropositives et 323 VIHséronégatives ont été criblés à le recherche d’ADN de VPH par PCR consensus au niveau du gène L1. Les échantillons positifs pour VPH-16 (n=74), -31 (n=74), -33 (n=37) et -35 (n=58) étaient analysés par PCR-séquençage pour la totalité du gène L1. Le nombre de nucléotides substitués pour L1 variait de 19 pour VPH-33 à 52 pour VPH-31. Le rapport du nombre de variantes sur le nombre d’isolats testés était plus élevé pour VPH-31 (56.4%, p=0.05) et VPH-35 (60.3%, p=0.04) comparativement à VPH-16 (40.5%), alors que ce ratio était inférieur pour VPH-33 mais sans différence statistiquement significative (24.3%, p=0.14). La distance entre les variantes était plus grande à l’intérieur des cinq boucles présumément exposées à la surface de la protéine L1 que dans la séquence à l’extérieur (p<0.01) Des variations synonymes étaient observées chez 1.7% (95% CI 1.1- 2.3) des nucléotides intra-boucles et 2.4% (95% CI 1.2-3.7) de ceux extra-boucles. Les variations non-synonymes étaient rencontrées pour 1.8% (95% CI 1.1-2.5) des nucléotides intra-boucles et 0.2% (95% CI 0-0.4) pour les nucléotides extra-boucles. Les ratios dN/dS étaient inférieurs à 1.0 pour les régions extra-boucles et encore davantage pour les régions intra-boucles. Ces résultats suggèrent que les séquences des régions hypervariables de L1 ont été sélectionnées positivement.

MiR-16, un nouveau régulateur du transporteur de glucose dépendant de l’insuline GLUT-4

Les microARNs sont des petits ARNs non codants d'environ 22 nucléotides qui régulent négativement la traduction de l'ARN messager cible (ARNm) et ont donc des fonctions cellulaires. Le microARN-16 (miR-16) est connu pour ses effets antiprolifératifs. Nous avons observé que l’expression de miR-16 est diminuée dans les cellules endothéliales humaines sénescentes et quiescentes en comparaison à des cellules prolifératives. Une analyse informatique des sites potentiels de liaison de miR-16 prévoit que GLUT-4, un transporteur du glucose insulinodépendant, pourrait être une cible potentielle du miR-16. Nous avons donc testé l'hypothèse que miR-16 régule négativement le métabolisme du glucose cellulaire. Dans des HUVEC, l'inhibition de miR-16 endogène avec des anti-miRNA oligonucléotides (AMO) augmente les niveaux protéiques de GLUT-4 de 1,7 ± 0,4 fois (p=0,0037 ; n=9). Dans des souris nourries avec un régime alimentaire normal ou riche en graisse et en sucre, l’expression de GLUT-4 dans le muscle squelettique a tendance à corréler négativement avec les niveaux de miR-16 (p=0,0998, r2=0,3866, n=4). Ces résultats suggèrent que miR-16 est un régulateur négatif de GLUT-4 et qu’il pourrait être impliqué dans la régulation du métabolisme cellulaire du glucose.

Sélection et polymorphisme chez des grenouilles mimétiques du Pérou (Dendrobatidae)

La diversification des signaux aposématiques dans un cadre de mimétisme müllérien est un phénomène intrigant. Alors que la théorie relative à l'aposématisme et au mimétisme suggère l'évolution vers un signal aposématique unique, d'impressionnantes variations peuvent être observées entre les populations, et cela à petite échelle spatiale. Il a été supposé que la variation spatiale des pressions de sélection engendrées par différents prédateurs puisse être à l'origine de ce phénomène. Afin de tester cette hypothèse, nous avons étudié la transition entre deux systèmes géographiques caractérisés par des patrons aposématiques distincts chez des grenouilles mimétiques et toxiques du nord du Pérou (Dendrobatidae) en combinant les outils de génétique des populations aux outils écologiques. Dans chacun de ces systèmes, Ranitomeya imitator vit en sympatrie avec R. ventrimaculata ou R. variabilis. Il s'agit du principal exemple empirique suggérant que dans un cadre de mimétisme müllérien, il n'y a pas convergence des signaux aposématiques des deux espèces, mais plutôt convergence unidirectionnelle où R. imitator, étant polymorphe, imite des espèces monomorphes avec lesquelles elle est sympatrique. Premièrement, les résultats réfutent les prémisses qui suggèrent que R. imitator converge vers le signal aposématique d’une autre espèce. La haute similarité génétique entre les espèces modèles suggère qu'elles ont divergé plus récemment que les populations de R. imitator ou qu'elles sont encore connectées par du flux génique. Ces résultats indiquent que ces espèces ont été identifiées à tort comme des espèces différentes. De fait, l'identification de l'espèce imitatrice basée sur la variabilité phénotypique est invalidée dans ce système puisque R. imitator et R. variabilis/ventrimaculata démontrent la même variabilité. Deuxièmement, nos résultats démontrent que la prédation varie spatialement, autant en intensité qu'en direction, créant ainsi un paysage hétérogène de pressions de sélection. Ainsi, de fortes pressions de prédation stabilisatrice permettent le maintien de l'organisation géographique de différents signaux aposématiques et expliquent l'uniformité de ces signaux ainsi que les relations mimétiques. Par contre, le relâchement temporaire des pressions de prédation permet l'apparition de nouveaux phénotypes aposématiques via les processus évolutifs neutres, conduisant à un haut polymorphisme au niveau de ces populations. L'interaction de ces modes sélectifs nous a permis de démontrer pour la première fois comment la théorie évolutive de Wright (shifting balance theory) permet la diversification adaptative dans un système naturel. Pour conclure, cette étude a permis de mettre en évidence à quel point les systèmes de mimétisme müllérien peuvent être dynamiques. L'alternance spatiale entre les processus évolutifs neutres et la sélection naturelle permet l'émergence de nouveaux phénotypes aposématiques à une échelle locale, ainsi que l'apparition d'une organisation géographique des signaux d'avertissement et des relations de mimétisme müllérien.